arXiv:1503.08460vl [cond-mat.mtrl-sci] 29 Mar 2015 


Graphene-like Dirac states and Quantum Spin Hall Insulators in the square-octagonal 

MX2 (M=Mo, W; X=S, Se, Te) Isomers 


Yan Sun^, Claudia Felser^, and Binghai Yan^’^ 

^ Max Planck Institute for Chemical Physics of Solids, 01187 Dresden, Germany and 
^ Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany 

(Dated: March 31, 2015) 

We studied the square-octagonal lattice of the transition metal dichalcogenide MX 2 (with M=Mo, 
W; Y=S, Se and Te), as an isomer of the normal hexagonal compound of MX 2 . By band structure 
calculations, we observe the graphene-like Dirac band structure in a rectangular lattice of MX 2 with 
nonsymmorphic space group symmetry. Two bands with van Hove singularity points cross each at 
the Fermi energy, leading to two Dirac cones that locates at opposite momenta. Spin-orbit coupling 
can open a nontrivial gap at these Dirac points and induce the quantum spin Hall (QSH) phase, the 
2D topological insulator. Here, square-octagonal MX 2 structures realize the interesting graphene 
physics, such as Dirac bands and QSH effect, in the transition metal dichalcogenides. 
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Since the discovery of graphene [TJ |2], research about 
two dimensional ( 2 D) materials has been widely explored 
in both theory and experiment. During the the last 
few years, transition metal dichalcogenides (TMD) MX2 
(with M=Mo, W and X=S, Se, Te) [3 [ 4 ] have attracted 
extensive attention. For example, monolayer of M0S2 
is a direct-gap semiconductor and regarded as massive 
Dirac systems in the honeycomb lattice with interesting 
valley physics [ 5 ]- Very recently an isomer structure of 
MX2 in the square-octagonal (so) lattice [TT] was found 
to exhibit gapless band structure with a Dirac cone a 
the zone center, and lattice distortion [12] was claimed 
to remove above Dirac cone and induce additional band 
crossing at the Fermi energy. However, the topological 
feature in the band structure was neglected due to the 
missing of spin-orbital coupling in calculations. In this 
work, we revisited the square-octahedral lattices of MX2 
isomers, and discovered their graphene-like Dirac band 
structures (see Fig. 1 ) and the 2 D topological insulator 
phase |6l Q • 

Inspired by the atomic structures of grain boundaries 
in normal hexagonal M0S2 iHin], a square-octagon (so) 
lattice o for M0S2 im was investigated in theory. There 
are four Mo atoms and eight S atoms in each primitive 
cell, as shown in Fig.l. And this lattice can be viewed 
as repeated square-octagon pairs in both x and y direc¬ 
tions. Here, the trigonal prismatic structure of M0S6 is 
slightly distorted compared to the hexagonal phase. Fur¬ 
ther, distortions from the square lattice to a rectangle 
lattice (distorted-so lattice) were found to optimize the 
strain in the 2 D structure and realize more stable struc¬ 
tures [12]. In this work, we adopted the similar lattice 
structures to TMD MX2 (with M=Mo and W; X= S, 
Se and Te) monolayers and found that the distorted-so- 
lattice is indeed the energetically favored structure for all 
compounds. The lattice distortion and total energy dif¬ 
ferences are summarized in Table I. We can see that for 
lattice distortion, characterized by the ratio of in-plane 


lattice parameters b/a, becomes stronger as the increas¬ 
ing of atomic radius of atom X for fixed metal atom M, 
due to the elongated M-X bond. 



FIG. I. (color online) Crystal lattice structure for MX 2 (M 
= Mo, W; X = S, Se, Te) in the (a) normal hexagonal lattice, 
(b) square-octagonal (so) lattice and (c) distorted-so lattice. 
The primitive unit cell is denoted by dashed lines. The band 
structures near the Fermi energy are illustrated in the 2D 
Brillouin zone for (d) so lattice, (e) distorted-so lattice. The 
3D plot of the distorted-so lattice is shown in (f), where two 
Dirac cones forms in the T — X line. The Dirac point and 
the van Hove singularity (VHS) point are indicated by empty 
and filled circles in (e), respectively. 

We mainly focus on their topological electronic prop¬ 
erties, after clarifying the lattice structures of MX2. 
Our calculations have been performed by using density 
functional theory (DFT) with projected augmented wave 
(PAW) method as implemented in the code of Vienna 
Ah initio Simulation Package (VASP) [131 Ej- The 
exchange-correlation energy are considered in the gener- 
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TABLE L The lattice constants and electronic properties of MX 2 (with M=Mo, W; A=S, Se and Te). Total energy differences 
between two types of lattices is defined as AE = Eso-iatuce — Edistorted-so-iatuce in one primitive cell. In the distorted-so- 
lattice, the linear band crossing at Fermi energy induces a Dirac cone, which locates on the line of T-X in the lattice momentum 
space. Its detailed location (kD) is defined as the the relative distance away from T point, as the schematic diagram given in 
Fig.1(b). The ab —initio calculations for lattice structure optimization, total energy differences AF, and location of Dirac cone 
are preformed without the inclusion of SOC, while for band gap Eg and Z 2 invariants the SOC was included. 



lattice constant (A) 

AE (meV) 

E, (meV) 

Z 2 invariant 

ko (f) 


so-lattice 

distorted-so-lattice 


so-lattice 

distorted-so-lattice 

so-lattice 

distorted-so-lattice 


M 0 S 2 

a=6.34 

a=6.30, 6=6.37 

0.68 

25 

12 

1 

1 

0.025 

MoSe2 

a=6.62 

a=6.56, 6=6.67 

1.11 

38 

23 

1 

1 

0.073 

MoTe2 

a=7.06 

a=6.72, 6=7.32 

17.87 

49 

19 

1 

1 

0.201 

WS 2 

a=6.36 

a=6.32, 6=6.42 

1.26 

111 

64 

1 

1 

0.072 

WSe2 

a=6.64 

a=6.30, 6=6.88 

29.50 

152 

20 

1 

1 

0.225 

WTe2 

a=7.11 

a=6.66, 6=7.37 

87.20 

213 

19 

1 

1 

0.245 


alized gradient approximation (GGA) level with Perdew- 
Burke-Ernzerhof (PBE) based density functional [ 15 ]. 
The energy cutoff was set to be 350 eV. The tight 
binding matrix elements were calculated by projection 
Bloch states onto maximally localized Wannier functions 
(MLWFs) [T7HT9] . using the VASP 2 WANNIER 90 inter¬ 
face [ 20 ] . 

Because all the compounds of MX2 share similar elec¬ 
tronic properties, in the following part we will take M0S2 
and WS2 as the examples for detailed analysis of their 
electronic structures. Band structures for two types of 
lattices are compared in Fig. 2 . Without the inclusion of 
SOG effect, band structure in the so-lattice presents dou¬ 
bly degenerated dz2 states at F point around the Fermi 
energy for both M0S2 and WS2, exhibiting as a semi- 
metallic state. As long as the SOC is taken into con¬ 
sideration, the degenerated dz2 bands split into two sin¬ 
gle states, with one locating at the top of valence band 
and the other up shifting to the second conduction band. 
Beside that, band anti-crossing between conduction and 
valence bands appears around F point, which implies the 
existence of band inversion. 

In order to make clear the topological electronic prop¬ 
erties, we analyzed the wave functions around the Fermi 
energy. As presented in Fig. 2 (a) and (d), the top of va¬ 
lence band and bottom of conduction bands at F point 
are mainly dominated by M — dz2 and M — dx2-y2 or¬ 
bitals, respectively, and these two states have opposite 
parities. Since the splitted two M — dz2 states are both 
with plus parities, the effect of SOC here is just opening 
the band gap but not changing topological band order. 
Because of the inversion symmetry in the so-lattice, we 
can directly achieve the topological number of Z2 invari¬ 
ant by the products of parities at time reversal invariant 
momentas (TRIMs) [T 6 |. The parity products for occu¬ 
pied states at the three independent TRIMs of r(0,0), 
X (( 0 . 5 , 0 ), ( 0 , 0 . 5 )) and S'( 0 . 5 , 0 . 5 ) are -, + and -h, re¬ 
spectively, which gives the Z2 invariant uq = 1 . There¬ 
fore, it directly conforms the existence of QSH insulator 
state in the S0-MX2. 



FIG. 2. (color online) Electronic band structures and density 
of states (DOS) for M 0 S 2 and WS 2 . Energy dispersion along 
high symmetry lines in 2D BZ for (a, b) M 0 S 2 with so and 
distorted-so lattices, and (d, e) WS 2 with so and distorted-so 
lattices, (c) and (f) are total DOS for distorted-so M 0 S 2 and 
WS 2 , respectively. Energy bands with (black solid lines) and 
without (red dot lines) are both included. DOS are calculated 
with the inclusion of SOC Local band structures around Fermi 
energy are also given in the inserted figures. VHS points in 
the DOS are remarked with red and blue dots. Fins and minus 
signs remarked in band structures are the parity eigenvalues. 
Band structures are calculated from ab — initio method, and 
DOS are calculated from MLWFs. Fermi energy is set to be 
zero. 

After the lattice distortion, the symmetry decreased 
from D2h to C2h^ and one can expect a dramatically 
change of the electronic properties. As shown in Fig. 2 (b) 
and (e), doubly degenerated M — dz2 states at F point 
split for both M0S2 and WS2, even without the inclu¬ 
sion of SOC. Meanwhile, a linear band crossing forms a 
massless Dirac cone at the Fermi level near F point on 
the line of V-X. The crossed two bands are mainly dom¬ 
inated by M — dz2 and M — dx2-y2 orbitals, respectively. 
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Besides, it is found that the location of the 2D Dirac 
cone kn strongly related to the distortion strength. As 
given in Table I, the Dirac cone shifts far away from T 
point as the increasing of lattice distortion. With the 
inclusion of SOC, as presented by the local band struc¬ 
tures in Fig. 2 (b) and (e), a general gap is opened with 
the breaking of the band crossing, which is just a typi¬ 
cal image for TIs. Since the distortion does not change 
the bulk band order, electronic structures are topologi¬ 
cal equivalent in so-lattice and distorted-so-lattice. For 
further confirmation, we also calculated the Z2 invari¬ 
ant. Though the distortion changes the lattice structure 
and atomic positions, inversion symmetry is preserved. 
So parity product at TRIMs is still effective for identi¬ 
fying the topological order [ 16 ]. Our calculations found 
that the Z2 invariant is 1 for any compound of MX2, 
as shown in Table I. Therefore, distorted-so-MX2 is still 
locating at topological non-trivial state. 

The SOC effect becoming stronger along with the in¬ 
creasing of atomic weight, and correspondingly, the SOC 
opened band gap should be also increased. As given in 
Table I, it is really the case in so-lattice. However, for the 
distorted-so-lattice, the lattice distortion is also becom¬ 
ing stronger as the increasing of atomic radius. Mean¬ 
while, the band crossing point kjj shifts far away from 
F point, as shown in Table I. Hence, the bulk band gap 
in the situation with distorted-so-lattice is decided by a 
competition between the strengths of SOC and lattice 
distortion. As presented in Table I, the competed result 
in this series of compounds gives the largest band gap of 
about 64 meV, appearance in WS2. 

Because of the non-trivial 2D bulk band order in MX2 , 
topological protected metallic edge state happens. For 
calculation of the edge state, we have constructed the slab 
model through MLWFs based tight binding method [ 13 ^ 
US]. In so-lattice, since the bands around Fermi level 
are almost consisted by the hybridized M — dz2 and 
M — dx2+y2 orbitals for all the compounds of MX2 , ML¬ 
WFs are derived from atomic dz2 and da,2-7/2-like orbitals. 
While for distorted-so-lattice, due to the difference be¬ 
tween lattice constants of a and 6, it is not accurate 
enough to describe the tight binding model by only in¬ 
cluding dz2 and (ia,2-y2-hke orbitals, in which da^^-like 
oribtal is also necessary. The tight binding parameters 
are determined from the MLWFs overlap matrix. The 
slab model was constructed in x direction for so-lattice 
due to the cubic symmetry, and both of x and y direc¬ 
tions were choose for distorted-so lattice, corresponding 
projected ID-BZ given in Fig.l (e). In order to elimi¬ 
nate the coupling between two edges, the widenesses of 
the slabs were up to 200 and 100 unit cells for M0S2 and 
WS2, respectively. 

From Fig .3 we can see that edge states exist for both 
so-lattice and distorted-so-lattice. While the details are 
depending on different compounds and edge termina¬ 
tions. For example, edge bands cut Fermi level three 


times for S0-WS2, whereas other cases just cut Fermi 
level once. As the differences of chemical potentials, some 
cases do not show Dirac point like edge states, but the 
non-trivial Z2 invariant guarantees the edge bands always 
cutting Fermi level odd times. 

(a) M0S2 (so lattice) (b) MoS^Cdistorted-so lattice) (c) MoS^Cdistorted-so lattice) 


(d) WS (so lattice) (0 WS2(distorted-so lattice) (f) WS2(distorted-so lattice) 

X r X 

FIG. 3. (color online) Tight binding edge band structures for 
(a) SO-M 0 S 2 , (b) distorted-so-MoS 2 , (c) so-WS 2 ,as well as (d) 
distorted-so-MoS 2 . Local band structures for M 0 S 2 around 
Fermi energy are showing inset. Red curves corresponding to 
edge states. Fermi energy is set to zero. 

As we have discussed above, the effect of SOC is just 
opening the band gap but not inducing band inversion. 
Therefore, the physics of band inversion can be under¬ 
stood without SOC. As in the situation of graphene, the 
band inversion in MX2 is due to the symmetry protected 
Van Hove singularity (VHS). In two dimensions, energy 
band with opposite energy dispersions along different di¬ 
rections can induce the saddle point at the transition 
momenta, which known as VHS. Moreover, if two bands 
with VHS cross each other, a band inversion will be there 
naturally. This is just case in the MX2. 

From the DFT band structures we have known that, 
without inclusion of SOC the so-lattice-MV2 presents as 
semi-metallic state with degenerated M — dz2 locating 
at F point. This degeneration can be regarded as an 
overlap of two saddle points belonging the two bands 
with VHSs, as the schematic showing in Fig. 1(d). After 
lattice distortion, the character of VHS becomes more 
obviously due to the symmetry breaking. As presented 
in Fig. 1 (e), the two neighboring bands shifts upwardly 
and downwardly, respectively, and overlapped two saddle 
points are separated by a band gap around F point. The 
features of VHS are also obviously in the local DFT band 
structures as given in Figs. 2 (a, b) and (d, e). Besides, 
according to the density of states (DOS) expression in 
lattice momentum space, g{E) oc f (where, g{E) 
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is the DOS at energy and i^j=x^y)^ DOS integrand 
at the saddle point diverges due to the extreme of energy 
dispersion. The diverged peaks can be found in the DOS 
for both M0S2 and WS2, as presented in Fig 2 .(c) and (f). 

Furthermore, as we have seen in the DFT band struc¬ 
tures, the lattice distortion also shifts the degenerated 
point away from F point. Similar to graphene, the new 
degenerated point exhibits as a massless Dirac cone, and 
the only difference is the location difference of the Dirac 
cone in lattice momentum space. In graphene, the Dirac 
cones paired locate at the high symmetry momentas of 
K and K' ^ which are connected by inversion and time 
reversal symmetry. While in distorted-so-MX2, Dirac 
cones paired locate at ko on the line of X-F-X, which 
are also connected by inversion and time reversal sym¬ 
metry, as shown in Fig. 4 (b) and (c). SOC is strong 
enough to open a considerable bang gap in distorted-so- 
MX2, different from the case of graphene with very weak 
SOC [2l]. 

In conclusion, by first principles calculations we have 
theoretically proposed a series of QSH insulators in the 
allotropes of monolayer TMD MX2 (with M=Mo, W; 
X=S, Se and Te ). The ground states in the allotropes 
show distorted-so-lattice. Similar to graphene, the band 
inversion is induced by symmetry protected VHS, and 
SOC just plays the role of opening band gap. How¬ 
ever, due to much stronger SOC effect in MX2, the in¬ 
verted band gaps are big enough for measurement. The 
band gap is about 12 meV in M0S2. And the largest 
band gap in this series QSH insulators is about 64 meV, 
which appears in in WS2. The QSH phase in MX2 
broaden the physical properties for DMTCs. As the so 
and distorted-so phase are derived from the grain bound¬ 
ary of h-MoS2 monolayers, it is possible to experimentally 
detect the QSH phase in M0S2 and the other TMD mate¬ 
rials through measuring the current or photo emission on 
their boundary. In addition, the VHS points indicate the 
existence of Lifshitz transition in the Fermi surface [ 22 ] 
and might promise interesting superconductivity in these 
materials. 

We are grateful to Z. Wang, C.-X. Liu, H. Su for fruit¬ 
ful discussion. 
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